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Abstract 

One of the open problems in nuclear structure is how to predict properties of finite nuclei from 
the knowledge of a bare nucleon-nucleon interaction of the meson-exchange type. We point out 
that a promising starting point consists in Dirac-Brueckner-Hartree-Fock (DBHF) calculations us- 
ing realistic nucleon-nucleon interactions like the Bonn potentials, which are able to reproduce 
satisfactorily the properties of symmetric nuclear matter without the need for 3-body forces, as 
is necessary in non-relativistic BHF calculations. However, the DBHF formalism is still too com- 
plicated to be used directly for finite nuclei. We argue that a possible route is to define effective 
Lagrangians with density-dependent nucleon-meson coupling vertices, which can be used in the 
Relativistic Hartree (or Relativistic Mean Field (RMF)) or preferrably in the Relativistic Hartree- 
Fock (RHF) approach. The density-dependence is matched to the nuclear matter DBHF results. 
We review the present status of nuclear matter DBHF calculations and discuss the various schemes 
to construct the self-energy, which lead to differences in the predictions. We also discuss how 
effective Lagrangians have been constructed and are used in actual calculations. We point out that 
completely consistent calculations in this scheme still have to be performed. 

PACS numbers: 21.30.Fc, 21.60. Jz, 21.65. Cd, 21.65.Mn 
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I. INTRODUCTION 

The Brueckner theory of nuclear many-body systems was introduced about 50 years ago[l| 
as a powerful tool to relate the bare nucleon-nucleon (N-N) interaction V to the effective 
in-medium interaction. Here and in the rest of this paper, we understand by bare interaction 
a meson-exchange type of interaction such as the Bonn potentials^] which are designed for 
describing two-body data up to a few hundred MeV. This effective interaction takes into 
account the effects of rescattering inside the nuclear medium and it is generally known as 
the Brueckner G-matrix. The G-matrix differs from the scattering T-matrix of an isolated 
N-N pair because of the presence of the surrounding nucleons which limits the available 
phase space. Thus, the G-matrix deviates from the bare interaction V, in particular by the 
fact that G behaves smoothly whereas V is strongly repulsive at short distances. Then, 
it becomes conceivable that a nuclear system can be described by an ensemble of nucleons 
moving in a mean field created self-consistently by the two-body effective interaction G while 
this picture would not hold in a simple Hartree-Fock (HF) description based directly on V. 
The Brueckner-Hartree-Fock (BHF) theory [3], is thus an essential link between the bare 
N-N interaction V and the properties of nuclear matter and atomic nuclei. 

n 

However, it was found already in early BHF calculations [5[ that the equation of state 
(EOS) calculated in BHF theory does not lead to the empirical saturation point of symmet- 
ric nuclear matter ( ~-16 MeV at kp ~1.36 fm _1 ). Indeed, it is well-known that all BHF 
predictions of the nuclear matter saturation point lie along the so-called Coester band (J 
which misses the empirical saturation point. It is only when contributions of 3-body forces 
are taken into account that the saturation point can be satisfactorily predicted [3, |8j. A 
large part of these 3-body terms were later understood as effects of the vacuum, i.e., vir- 
tual nucleon-antinucleon excitations, leading to the so-called Z-graphs in the diagrammatic 
expansion 

an 

It is then reasonable to include effects of relativity from the beginning. The relativistic 
version of the BHF theory is known as the Dirac-Brueckner-Hartree-Fock (DBHF) approach 
since the early work of Refs. ll|, |l2[ and the fully self-consistent nuclear matter calculations 



of Ter Haar and Malfliet 



13] . The starting point of the DBHF theory is a realistic bare N-N 



interaction formulated in terms of meson exchanges, the parameters of which are adjusted to 
reproduce N-N scattering data in the T-matrix approximation. This One-Boson-Exchange 
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(OBE) interaction is used in a relativistically covariant framework by solving the Bethe- 
Salpeter (BS) equation, or more precisely a Thompson equation obtained by a covariant 
three-dimensional reduction of the BS equation |l3j. In the relativistic approach the nu- 
clear mean field is obtained by a large cancellation between scalar and vector potentials. 
The scalar potential saturates at high densities, thus providing an additional saturation 
mechanism. It is then found that the DBHF approach based on meson-exchange potentials 
predicts a nuclear matter saturation point away from the Coester band without the need 
for a 3-body force and in relatively good agreement with the empirical point HQ. This 
advantage of DBHF over the non-relativistic BHF makes the applications to finite nuclei 
practicable in principle and, in fact, more simple. 

The Brueckner theory is technically difficult, which is related to the double self- 
consistency of the problem: the determination of the effective in-medium interaction G 
requires the knowledge of the self-consistent single-particle mean field, but this mean field 
itself depends on G. This double self-consistency can be handled in simple cases like homoge- 
neous infinite nuclear matter but it becomes very cumbersome and difficult to treat in finite 
nuclei. There are some BHF calculations of finite nuclei, e.g., Ref. 15[, but not including 
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3-body forces, and only very rare calculations in the relativistic DBHF version 

Then the question arises whether one can effectively perform self-consistent mean field 
calculations for finite nuclei based on the knowledge of the mean field properties of infinite 
nuclear matter deduced from BHF or DBHF calculations. The idea is to build an effective 
Hamiltonian that can reproduce the results of the non-relativistic BHF (plus corrections), or 
the relativistic DBHF in infinite nuclear matter in simpler many-body schemes, namely the 
non-relativistic Hartree-Fock (HF) or the relativistic Dirac-Hartree-Fock (DHF) approxima- 
tions, and then solve the HF or DHF e quat ions for finite nuclei. This issue was first taken 
up in the 1960's and the work of Negele 17j showed that it was indeed possible to describe 
finite nuclei reasonably well using a G-matrix coming from infinite matter calculations, at 
the cost of some readjustments which were attributed to the neglect of 3-body effects and 
higher order terms. 

In the light of the success of nuclear matter DBHF calculations it then seems most 
promising to look for effective Lagrangians which are written in terms of OBE interactions 
with density-dependent meson-nucleon coupling vertices. The nuclear matter DBHF results 
would have to be reproduced in DHF approximation by these Lagrangians. Then, these 
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adjusted Lagrangians would be used for describing finite nuclei in the DHF approximation. 
So far, this program has never been carried out, and this constitutes the open question which 
is addressed in this paper. 

This scheme is in the spirit of a density functional approach, since it leads to an energy 
functional in terms of the baryonic density alone, which contains the effect of exchange and 
correlations. The application of such a Lagrangian to finite nuclei implies a sort of Local 
Density Approximation (LDA). This has been justified in noticing that the range of the 
effective interactions due to meson exchanges are smaller that the scale of spatial variations 
in a nucleus [13]. This assumption is weakest for 7r-meson exchange, but this field has special 
properties anyway, because of its pseudo-scalar (PS) nature. It should also be understood, 
that these meson exchanges do not necessarily represent the exchange of physical in-medium 
mesons, but they rather represent the different Lorentz components of the decomposition 
of the DBHF self-energy, such as scalar /vector or isoscalar/isovector fields. One should 
therefore rather speak of meson-like fields. In this spirit the masses of these various mesons 
do not necessarily have to be the masses of the physical, even less the free, mesons, even 
though this has been the usual assumption. 

The density dependence of the effective Lagrangian (i.e. of the coupling vertices) is 
determined by comparison of the relativistic Hartree (also called the Relativistic Mean Field 
(RMF)) or the relativistic Hartree- Fock (also called Dirac HF (DHF)) approaches with the 
full DBHF results. The question is then to what physical quantities of DBHF they are 
matched. This question arises because in practice the DBHF equations are often solved in 
a three-dimensional Thompson reduction, also neglecting the negative energy states (no-sea 
approximation). Then, the Lorentz structure of the self-energy is either lost or ambiguous. 
This will be discussed in the next section. 

The infinite matter DBHF calculations, on which this scheme is based, have limitations 
at very low and very high densities. At very low densities they become instable because 



of the deuteron po 
of nuclear matter 



e, or more generally because of cluster effects and non-inhomogeneities 



18J. At densities several times the value po of saturation density the 



description in terms of nucleons and mesons may not be valid any more, and more nucleon 
and meson states (also strangeness degrees of freedom), and eventually subnuclear degrees 
of freedom have to be included. This is the regime of more or less energetic heavy ion 
collisions, where these questions are also of immediate relevance. They are also of relevance 
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in the structure of neutron stars and core-collapse supernovae, where the question of the 
density dependence of the symmetry energy enters in particular. However, in the regime of 
nuclear structure, i.e., in the range of about 0.1p to (1 — 2)p the approach sketched above 
may be valid. It is this approach which we want to review in this article. 

The structure of this article is the following. In Sect. II we shall review the status of 
DBHF calculations for infinite nuclear matter, in particular the problem of constructing the 
various Lorentz components of the self-energy, also discussing asymmetric nuclear matter. In 
Sect. Ill we discuss the question of constructing effective (density- dependent) Lagrangians 
from the DBHF solutions in nuclear matter. We also make sime comments on applications 
of such effective Lagrangians to finite nucleus calculations. Concluding remarks are made 
in Sect. IV. We hope it will be apparent, that on one hand the scheme of effective DBHF 
Lagrangians is an important and viable route from the realistic N-N interactions to realistic 
calculations of finite nuclei. Compared to the highly successful phenomenological RMF and 
RHF calculations this approach is essentially parameter-free, in the sense that there are no 
parameters adjusted to nuclear matter or finite nuclei. On the other hand, we will see that a 
completely satisfying solution to this problem is still to be found. 



II. STATUS OF DBHF CALCULATIONS IN INFINITE MATTER 



The earliest study of infinite nuclear matter based on the self-consistent framework of 
DBHF theory was performed by Erkelenzj^j. As his work also explored the description of 
the N-N scattering data using a one-boson exchange (OBE) interaction, he compared the 
results obtained using several early versions of the OBE without attempting to reproduce the 



Hi, 



expected saturation properties of nuclear matter. With the success of Dirac phenomeno 
in the description of proton-nucleus scattering, the problem was taken up again in Ref. 
where reasonable agreement with the expected nuclear matter properties were obtained 
using the OBE N-N interactions available at the time. In the pioneering work of Ter Haar 
and Malfliet 13|, a one-boson exchange interaction reproducing N-N and N-A data in the 



energy range up to 1 GeV was used in DBHF calculations of symmetric nuclear matter 
and the effects of the N-A degree of freedom were also studied. In the meantime DBHF 
calculations have reached a high level of sophistication. They confirm that, starting from 
high precision N-N interactions, one may obtain good predictions for nuclear matter. This 
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also gives confidence that this approach can be used to extrapolate to extreme states of 
nuclear matter at high density and temperature or high asymmetry. Such predictions are 
important for investigations in exotic nuclei, heavy ion collisions, and astrophysics. 

However, DBHF calculations do not yield identical results, even when the same N-N 
interactions are used. This has its origin in the difficulty mentioned earlier of the requirement 
of double self-consistency of both the G-matrix and the single-particle mean field in solving 
the Bethe-Salpeter equation. To simplify its solution, the reduced Bethe-Salpeter equation is 
usually projected on the positive-energy Dirac spinors, which yields an equation very similar 
to the non-relativistic one and permits a similar solution. However, the Lorentz structure of 
the G-matrix is lost in such a reduction and must be reconstructed in some manner. Various 
approximation schemes have been used to deal with this problem. This not only affects the 
results of the DBHF calculations itself, but also the construction of effective Lagrangians to 
be used for finite nuclei calculations, as discussed in the next section. 

From parity and time reversal symmetry the most general form of the self-energy in 
nuclear matter has the form 

E(Jfe) = E s (fc) -7°S°(fc) +7- k £"(fc) , (1) 

with the scalar, time and space-like vector components of the self-energy. In asymmetric 
nuclear matter these are different for protons and neutrons and the appropriate isoscalar and 
isovector components have to be determined in addition. In the Hartree and also in the HF 
approximation these follow in a unique way from the exchange of the corresponding meson 
fields. In DBHF, however, the exchange and correlation diagrams mix the contribution 
of the various mesons to the components of the self-energy. Thus, 7r-meson exchange, for 
example, contributes to the self-energy even though there is no pseudo-scalar (PS) self-energy 
in Eq.([T]). The above remark also explains why in the Introduction we made the comment 
that one should talk of meson-like DBHF fields and not of the exchange of physical mesons. 
Thus, the Lorentz structure has to be reconstructed, giving the correct contribution of the 
various meson exchanges. 

The simplest method to determine the above components of the self-energy was intro- 
duced by Brockmann and Machleidt jl4, [jjj and consists in analyzing the momentum de- 
pendence of the single-particle energies, 

e(k) = Vk 2 + M* 2 + S , (2) 



where the effective mass is M* = M + S s . It was later found 20|, |2l| that, while this 
method gives reasonable results for symmetric nuclear matter, it leads to incorrect results 
in asymmetric matter, even giving the wrong sign for the isospin dependence of the effective 
interaction. The alternative approach is to project the G-matrix elements on a complete 
set of five - in asymmetric matter six - Lorentz invariants 22j. However, the choice of these 
invariants is not unique. Different choices give identical matrix elements for the Hartree self- 
energies for on-shell positive energy states, but the exchange and correlation contributions, 
as well as the coupling to negative energy states, can be very different leading to different 
self-energies. Thus, various projection schemes have been developed which differ mainly in 
the effect of the PS meson exchange. 

One possible approach is to use the PS operator as one of the invariants. The pseudo- 
vector (PV) pion-nucleon coupling is defined so as to yield the same matrix elements between 
positive energy states as the PS coupling. A PS coupling, however, furnishes very strong 
coupling between positive and negative energy states and it is unsuitable for use in nuclear 
matter calculations. As the covariant reconstruction cannot distinguish between the PS and 
PV forms, it can furnish unexpectedly large and undesired PS contributions. It also leads 



to a very strong momentum dependence of the self-energies, see e.g., Ref.[23|. Therefore, 
the preferred choice has been to use a pseudovector PV invariant. However, this is also not 
completely satisfactory, because the p-exchange for instance is not treated correctly then 23]. 
A compromise, which has proven to be most successful, is a decomposition of the DBHF 
G-matix of the form 

G = V + AG , (3) 
where V is the starting OBE interaction. This approach was introduced by the group of 



Miither 



20 



24j, and in a slig htly different form as the 'optimal T-matrix representation' 



by the group of Fuchs 



23 



25 



26j]. Since the Lorentz structure can be derived explicitly in 



the HF approximation for the OBE interaction, the first part of V + AG is treated exactly. 
Only the difference AG is treated by a projection method. This method seems to yield the 
most reliable results at present. In this method AG, which represents mainly the correlation 
effects, has been parametrized in terms of effective meson exchanges of small or zero ranges 
with density dependent coupling vertices. In a related approach Boersma et al. 27|, |28j have 
parametrized the G-matrix in terms of Yukawa function with fitted masses and coupling 
strengths including a form factor. These effective interactions have then been used in finite 
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nucleus calculations directly (see next section). 

The PV projection approaches with or without splitting off the bare OBE interaction 
usually agree rather well for symmetric nuclear matter around saturation density. They give 
differences, however, for asymmetric or neutron matter, and also for higher densities. Thus, 
these different choices are of importance for the structure of exotic nuclei with neutron-rich 
low-density skins and also for heavy ion collisions and the structure and mechanisms of 
astrophysical objects. 

Still another method consists in performing the projection on the full space of positive 
and negative energy Dirac states, as proposed e.g. by de Jong and Lenske 29J and also 
in a slightly different framework by Huber at al. 30(. There would then be no ambiguity. 



However, the authors of Ref. 



29j point out that, in order for this scheme to be satisfactory, 



the N-N T- matrix would have to be determined in the same scheme, which is not usually 
the case. Thus, this method also needs further development. 

It is well known that the DBHF is not a thermodynamically consistent approximation 
3l| . Several attempts have been made to go beyond it by developing a Dirac T- matrix 
approximation to nuclear matter 32|, |33[. These yield results in agreement with previous 
ones near saturation but tend to predict a stiffer equation of state at higher densities. A 
calculation explicitly including the negative energy states has also been performed 29J and 
it also yields a stiffer equation of state at high densities. However, these calculations require 
approximations unnecessary for the DBHF, such as the quasi-particle approximation, and 
should be improved before conclusions are made about their high density behavior. 

Other approximations and difficulties are relevant to the DBHF. Schiller, Miither and 
Czerski 34( have shown that binding energies near saturation in Brueckner calculations can 



vary by about 10% due to angle averaging of the Pauli blocking factor. Clustering creates 
serious problems for obtaining convergent solutions at low nuclear matter densities and, even 
where convergent, furnishes binding energies that cannot be described by a Hartree-Fock 
type calculation 35| . At sufficiently high densities, one can also question the validity of a 
OBE approximation to the interaction, as discussed in the Introduction. Yet, in spite of the 
many open questions, the existing nuclear matter DBHF calculations indicate that one can 
describe the bulk properties of infinite matter in terms of realistic OBE interactions, at least 
at values of matter density relevant for nuclear structure calculations, which are neither too 
low nor too high. 
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III. DBHF-BASED EFFECTIVE LAGRANGIANS 



The next step is to determine effective Lagrangians which can be used in a self-consistent 
mean field approach like RMF or RHF and which reproduces the DBHF results in homo- 
geneous matter. The effective Lagrangians is described as a system of nucleons of Mass M 
and effective mesons interacting via meson-nucleon couplings which are determined by the 
coupling strengths gi in the following form 

£ = $ (17^ — M) ipi + C° mes (a, u, p, 1 7i, A) + C mt (4) 

£mt = 9ai><^> - 9<Jh,kU*i> - g P 4>i^ ■ rip + g s ^S- rip (5) 

/ _ 1 

H ^ ^757^7? ■f^-e^7Mo( 1 + r 3)^V' , 

m n 2 

with the fields of the nucleon (ip), two isoscalar meson fields (a and u), three isovector 
fields (7?, p and 5) and the photon (A). The free Lagrangian £^ es (cr, u, p, 5, tt, A) of the 
mesons and the photon has the usual form with masses rrij of the mesons, which are often 
assumed to be the masses of the corresponding free mesons, even though this does not 
have to be the case considering the effectiveness of the meson fields. In many models the 
5-degree of freedom is omitted and the isovector properties rely only on the p meson if the 
RMF (Hartree) framework is adopted. In this case the pion is also omitted since it cannot 
contribute to the Hartree energy due to parity conservation. On the other hand, in RHF 
all mesons contribute to isovector properties, even isoscalar mesons through their exchange 
contributions to self-energies. 

The stationarity condition of the action integral J d A xC(x) under variations of the physical 
fields <p (<f) = if), a, u, p, ft, 5, and A) leads to the Euler-Lagrange equations, from which one 



can derive the equations of motion for the meson, photon and nucleon fields (lQj. The 
meson and photon fields obey inhomogeneous Klein-Gordon equations and a Proca equation 
with source terms, respectively, whereas the nucleon field obeys a Dirac equation coupled 
minimally to the meson fields. In homogeneous matter they have simple analytic forms 
which make the solving of the self-consistent problem much easier than in finite nuclei. It is 
customary to make the simplifying assumption of neglecting the time component of the four- 
momenta carried by the mesons, i.e., the meson fields are assumed to be time independent. 
This amounts to neglect the retardation effects. The nuclear energies involved are small 
compared to the masses of the exchanged mesons, so that this approximation should be 
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valid for the a-, u-, 5- and p- induced interactions, and also, to a lesser extent, for the pion. 

The meson-nucleon couplings must be determined by reproducing as closely as possible 
the DBHF results of infinite matter in some range of values of baryonic density and neutron- 
to-proton ratio. Since the DBHF self-energies are density and momentum dependent, the 
coupling vertices necessarily have to be density dependent, (we will discuss the momen- 
tum dependence below), i.e., one has g%{pB) for the different meson species i, where ps is 
the baryon density. Of course, the various approximations schemes for the self-energies of 
DBHF in nuclear matter, which were discussed in sect. II, will lead to different results also 
for the effective Lagrangians. As discussed there, the method of extracting these self-energies 
from their momentum dependence [l^] was shown to give incorrect asymmetry properties 20]. 
Thus, we will concentrate on approaches which are based on matching the Lorentz compo- 
nents of the self-energy. Finally, there are also the approaches based on parametrizations of 
the DBHF G-matrix directly. 

Two scenarios have been applied with respect to the approximation in which the DBHF 
self-energies are matched: either in the Hartree or the HF approximation of nuclear matter. 
In the former case this has been called the density-dependent relativistic Hartree (or mean 
field) approximation (DDRH or DDRMF), in the second the DDRHF approximation. In 
DDRH the result is simple: the density- dependent coupling vertices are given via the ratio 
of the DBHF self-energies divided by the respective Lorentz components of the density, e.g., 

9a_\ 2 = 1 K( k F P , k Fn ) + Z s p (k Fp , k Fn ) 
mj ~ 2 p° n + p; 

and correspondingly for the other coupling coefficients (those for the 7r-coupling are given 



in ref.[36j). In Eq.([n]) we have explicitly indicated that the coupling vertices depend on 
the neutron and proton densities separately (resp. the isoscalar and isovector densities). 
In practice it was often found, that the dependence on the isovector density is weak and 
this dependence has been neglected (see, e.g., ref-jsil). In the Hartree- Fock approximation 
(DDRHF) the expressions which have to be matched are slightly more complicated, and are 



given also in ref. 36]. 



The DBHF self-energies are explicitly momentum dependent due to the contributions 
of the exchange and correlation terms (the exchange contribution follows roughly a 1/k 
dependence). In the DDRH approach, however, the RH self-energies in nuclear matter 
are momentum independent, and in the fitting procedure some approximation has to be 

11 



made. The approach followed, was either to match the self-energies at the Fermi surface 



(e.g. 



381 ] ) or to average over the Fermi sphere 36|. Another approach has been to perform a 



first order correction involving the derivative of the self-energies with respect to momentum 
at the Fermi surface {stJ. A covariant way to parametrize the momentum dependence in 
phenomenological RMF models has been suggested by Typel 18[ , but has not yet been used 
to parametrize DBHF self-energies. 

In the DDRHF approach the RHF expressions in nuclear matter are density dependent 
(roughly as 1/k) and thus the matching to the DBHF self-energies can be done in a more 
consistent way. In principle the DBHF self-energies also carry a momentum dependence 
from the correlation terms, but the matching to the 1/k dependence usually is quite accu- 



rate, see e.g. ref.[36|. Thus, the DDRHF approach should be preferred and should give a 
better connection between nuclear matter and finite nuclei. Fock exchange terms can be 
quite important in finite nuclei, as they can modify the shell structure or avoid unphysical 



shell gaps |c 



. RHF calculations in finite nuclei are more involved than RH ones but the 



additional technical difficulties can be solved 101 ] . 

Of course, the parametrization of the coupling vertices and the approximation in the 
calculation of finite nuclei has to be consistent. A first check has been that the equation 
of state, i.e., the energy per particle as a function of density, of DBHF is reproduced with 
the density- dependent parametrization in the same approximation (DDRH or DDRHF). A 
parametrization of DBHF in the Hartree approximation should be used in DDRH (DDRMF) 
calculations of finite nuclei, and the one obtained in RHF should be used with DDRHF 
calculations. 



Finally the G-matrix of DBHF calculations has been parametrized in terms o 



change of effective mesons by Schiller et al. 



24j or by Boersma and Malfliet 



27 



the ex- 



28 1 . These 



effective interactions have then been used directly in a finite nucleus calculation |40|, |4l| or 
used to reconstruct a self-energy which is used in turn to build a density-dependent effec- 
tive Lagrangian 42]. In this last procedure, however, the RHF method should be used for 
consistency. 

A Lagrangian as in Eq.([6]), in which the coupling strengths depend explicitly on the 
density p# violates the covariance of the model. This has been treated by writing the 
density dependence in terms of scalar invariants built from nucleon field operators related 



to the density 



431 ] . The form of this operator is not a priori given, and different choices 
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have been investigated, e.g., the scalar or the square of the vector density operator. The 
dependence on field operators necessarily leads to additional terms in the Euler-Lagrange 
equations of motion, so-called rearrangement terms, which depend on the type of density 



dependence assumed. In the investigations of Fuchs et al.[43j] it was found that a vector- 
density dependence gave better results for finite nuclei starting from the same effective DBHF 
self-energies, and such a dependence is now generally employed. In any case rearrangement 



contributions should always be included. It was shown in Ref.|43| that they are necessary 
for the thermodynamic consistency of the effective theory. The effect of the rearrangement 
terms in phenomenological density-dependent models has recently been studied in detail in 



Refs. 



44 



45|. 



Attempts have been made to calculate the properties of finite nuclei with the above men 
tioned effective Lagrangians 



or in the DDRHF approaches 



eluded in Refs 



hese studies have been done in the DDRMF 



21 



37 



38 



42 



28 



36 



43| 



37 



42 



4l|. Consistent rearrangement contributions were in- 



43[, but not in Refs. 



2N 



36 



38 



41| . From the arguments presented 



in this article a complete calculation should be done in DDRHF including rerrangement 
terms and matchin g to DBHF self-energies, which are obtained from a "subtracted" G- 



Matrix like in Refs. 



24 



26] . Looking through the calculations existing at present such a 



calculation does not yet exist, even though all the ingredients seem to be present. 

It is not the place here, to compare in detail the different results of the finite nucleus 
calculations based on DBHF calculations. Today, the focus is especially on the isovector 
properties and the extrapolation to very neutron-rich nuclei. Generally it can be said that 
the quality is comparable to the results of the phenomenological DDRMF or DDRHF cal- 
culations, where the density dependent parameters are adjusted directly to data of nuclear 
matter and finite nuclei with a number of (global) parameters ranging from 7 to 10, which 
naturally lead to a greater precision with respect to the data. Recently, the calculation 
scheme for a complete relativistic Hartree-Fock-Bogoliubov (RHFB) description of finite nu- 
clei has been proposed |46|, |47[ and the results are at an excellent quantitative level using a 
phenomenological Lagrangian. 
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IV. CONCLUSION 



In this work we have addressed the question of how to link the bare N-N interaction V 
to the properties of infinite and finite nuclear systems. This issue is of utmost importance 
because of the need to have a theory capable of predicting nuclear properties in regions 
of the nuclear chart not yet accessible to experiment. The self-consistent mean field is a 
good approach, but its predictive power cannot be extrapolated very far if one relies only 
on phenomenological adjustments of effective interactions or energy density functionals, as 
it is done in non-relativistic approaches. There are also attempts in non-relativistic mean 
field studies to use the bare N-N interaction limited to some finite part of momentum space, 
for instance with the VJ ow -k interactions. Here, we have discussed the possibility to deduce 
the effective N-N interaction via Brueckner theory. This could be done in a non-relativistic 
approach or a relativistic one. The main advantage of the relativistic Brueckner theory, or 
DBHF, is that the nuclear matter equation of state is satisfactorily obtained without adding 
explicitly the effects of 3-body forces as it is the case for the non-relativistic Brueckner 
theory. Thus, the DBHF equation of state and self-energies are good starting points for 
deducing effective Lagrangians to be used in relativistic self-consistent studies. 

With respect to the present state of the subject, there are not yet completely satisfactory 
RHF calculations of finite nuclei based on effective Lagrangians deduced from DBHF nuclear 
matter studies. The main difficulty was identified as coming from the ambiguity in the choice 
of the covariant operators on which the effective G-matrix has to be expanded. A way to get 
around has been proposed by Schiller and Miither 24| who suggest to split G into V + AG, 
the idea being that V is the dominant part of G and thus, a large part of the ambiguity can 
be eliminated. Indeed, there have some attempts to build G in this way 42| but the use of 



such parametrized, density-dependent G-matrix in RHF calculations of finite nuclei is still 
to be done. 

With the recent progress in RHF and RHFB calculations using density- dependent meson- 
nucleon couplings, the next step to establish a link between N-N interactions like the Bonn 
potentials and the properties of finite nuclei is at hand. This is one of the open problems 
in nuclear structure which need to be solved, and the route through DBHF and density- 
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dependent RHF seems a promising option. 
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